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Abstract 

We study Heisenberg model of classical spins lying on the toroidal support, whose 
internal and external radii are r and R, respectively. The isotropic regime is characterized 
by a fractional soliton solution. Whenever the torus size is very large, i? — > oo, its charge 
equals unity and the soliton effectively lies on an infinite cylinder. However, for R = the 
spherical geometry is recovered and we obtain that configuration and energy of a soliton 
lying on a sphere. Vortex-like configurations are also supported: in a ring torus (R > r) 
such excitations present no core where energy could blow up. At the limit R ^ oo we are 
effectively describing it on an infinite cylinder, where the spins appear to be practically 
parallel to each other, yielding no net energy. On the other hand, in a horn torus {R = r) 
a singular core takes place, while for i? < r (spindle torus) two such singularities appear. 
If R is further diminished until vanish we recover vortex configuration on a sphere. 

1 Introduction and motivation 

Geometrical and topological concepts and tools are important in many branches of natu- 
ral sciences, namely in Physics. For instance, the idea of symmetry, intimately associated to 
geometry, is a keystone for studying a number of fundamental properties of several physical 
systems, e.g., Noether theorem asserts a conserved quantity to each continuous symmetry of the 
associated action. Topology, in turn, is crucial for classifying and for giving stability to certain 
excitations, like solitons (extending objects having finite energy) and vortices (presenting non- 
vanishing vorticity around a given singular point or a topological obstruction). These and others 
have been observed in a number of systems, like superconductors, superfluids, and magnetic ma- 
terials, namely vortex-like magnetization has been directly observed in nanomagnets [Ij. On 
the other hand, the observed vortex-pair dissociation is the mechanism behind the topological 
phase transition [2]. Another kind of topological objects emerges from Euclidian non-Abelian 
pure-field models in (3+1) dimensions. They are called instantons, since they represent points 
in the Euclidian space-time. In this case, the relevant homotopy group is that associated to the 
mapping of the field space to the 4-dimensional Euclidian one, so that TCsiS^ — >■ S^) = Z, for 
example, if the internal space is the QCD gauge group, SU{3). Such excitations play important 
role in this framework once the nonconservation of the flavor singlet axial current is proportional 
to the instanton charge [3]. 
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In turn, the toroidal geometry has recently received a considerable attention, for exam- 
ple, as tight traps for Bose-Einstein condensates, inside which the atoms develop quasi-one- 
dimensional confined dynamics subject to periodic boundary conditions [5]. Ring-shaped carbon 
nanotubes have been shown to provide quasi-zero-dimensional systems, whenever the rings are 
very small[5l E]. In addition, their topology also make them possible field-effect transistors 
for technological applications [7j, as well as, the suitable control of their stable magnetization 
(e.g., those with vortex-like profile) may be useful for applications in magneto-electronic de- 
vices, namely if an array of magnetic nanorings is concerned [8]. From a more fundamental 
point of view, it has been verified that whenever the Ising model is defined on a donut-shaped 
lattice, then instead of an unique, there appear two critical temperatures [9]. Even in Biology, 
the toroidal shape plays important roles: for instance, it has been observed that a large number 
of proteins involved in DNA metabolism adopts a ring-like shape, even though these proteins 
have quite distinct and unrelated functions in this mechanism. Why this geometry/topology 
is so abundant in this process and the reason why life evolution has selected it amongst many 
others remains a puzzle in the specialized literature [lOj. 

Geometrically, a torus of genus-1 (one central hole) is a compact surface whose (Gaussian) 
curvature smoothly varies from —l/r{R — r) to +l/r{R + r), along its polar angle (see details 
below), so interpolating between the pseudospherical and spherical curvatures whenever R> r. 
From the topological point of view, the simplest genus-1 torus is obtained by the topological 
product of two circles, = A [H]. [In general, a genus-n torus is got in a similar way, for 
example, that of genus-2 is given by = A T^]. This confers it a non-simply connectivity, 
say, closed loops holding on its central hole or the polar circumference cannot be shrunk to a 
point (other non-simply connected topologies are the ordinary circle and a disk with a hole, an 
annulus) . 

For studying classical spin-like textures, of solitonic and vortex kinds, on a torus, we first 
write down the continuum version of the Heisenberg (exchange) Hamiltonian in this geometry 
(Section 2). Later, the non-linear Euler-Lagrange equations are obtained and some particular 
cases are explicitly considered for describing the desired excitations. In addition, their profiles, 
energies and other basic properties are discussed and compared to their counterparts from other 
surfaces. These tasks are performed in Section 3, where solitonic-like solutions in the isotropic 
regime are described and, in Section 4, which is dedicated to vortex-hke configurations, stud- 
ied within the Planar Rotator Model (equivalently, the XY model, once are dealing only with 
static properties). Finally, we close our paper by pointing out our Conclusions and Prospects 
for forthcoming investigation. 



2 The continuum Heisenberg model on the torus 

The anisotropic exchange Heisenberg model, for nearest-neighbor interacting spins on a two- 
dimensional lattice, is given by the Hamiltonian below: 

Hlatt = -J' = E ("^^"^i + ^^^J + (1 + ^'^^^^J^ 

<i,j> <ij> 

where J' denotes the coupling between neighboring spins, and according to J' < or J' > 0, the 
Hamiltonian describes a ferro or anti ferromagnetic system, respectively. Si = (Sf , Sf , S-) is the 
spin operator at site i, the parameter A accounts for the anisotropy interaction amongst spins: 
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R 

Figure 1: Shape of an ordinary ring torus embedded in a three-dimensional space (Left). Its cross 
section and the variables used throughout this work (Right). Azimuth-like angle, (p, runs along the 
torus tube, z-constant planes, and is not shown in the cut. 



for A > 0, spins tend to align along the (internal) Z-axis (easy-axis regime); forA = 0, one gets 
the isotropic case; for — 1 < A < 0, we have the easy-plane regime, while the A = —1 case, yields 
to the so called XY model (or the Planar Rotator Model, PRM, if we focus on 2-component 
spin, imposing = 0, so that Sprm = {S^, Sy)). 

In the continuum approach of spatial and spin variables, valid at sufficiently large wavelength 
and low temperature, the Hamiltonian ([T]) may be written as (J = J'/2) [T2[ ITSl ITi]: 



= jj^j{l + 6a3\){DS''fdQ (2) 

where f2 is the surface with curvilinear coordinates rji and 772, so that d^l = \/\g\driidri2, Sas 
is the Kronecker symbol, D is the covariant derivative, \/\g\ = >y\det[gij]\, gij and hat are 
the elements of surface and spin space metrics, respectively (as usual, gijg^^ = 5^). Now, 
S = {Sx,Sy,Sz) = (sin 9 COS sin 9 sin $, COS 9) is the classical spin vector field valued on a 
unity sphere (internal space), so that 9 = 9(771,772) and $ = $(771,772) . With this Carte- 
sian parametrization for S we have hab = Sab- The Hamiltonian ([2]) may be also viewed as an 
anisotropic non-linear a model (NLcxM) , lying on an arbitrary two dimensional geometry, so that 
our considerations could have some relevance to other branches like hydrodynamic, superfiuidity 
and superconductivity. 

Our interest is to study the model above on the torus geometry which is a smooth surface 
with varying curvature. The simplest torus is thought as a surface having genus one (a single 
central hole), and whenever embedded in three dimensional space it shapes like a donut (see 
Fig. [1]). The standard tori are classified in three distinct types concerning the relations between 
their internal, r, and external, R, radii. For i? > r, we have a ring torus (donut shape, as shown 
in Fig. [1]), if i? = r a horn torus is obtained, while for R < r a spindle torus is described (the 
latter ones are depicted in Fig. [2]). 

Any ordinary torus may be parametrized in a simple way by some distinct coordinate sys- 
tems, Uke the Cartesian [15] and the peripolar [16] (6*, 0), so that: 



(^R_,/^^:^y + z2=r\ (3) 
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with the parametric equations 



X 



y = {R + r sin 6) sin i 



r cos 9 



(4) 



{R + r sin 9) cos 

where R and r are the rotating (external) and axial (internal) radii, respectively (see Fig. [T]). 
In peripolar variables {9, (p) the metric elements read: 

(5) 



gee 



{R + r sin i 



and 



^ 9e<p = 0, 



from what there follows the Gaussian curvature: 



K 



sin 9 



r{R + r sin i 



(6) 



Note that K varies from -l/r{R-r), at 9 = 37r/2, to +l/r(i? + r), at 6* = 7r/2. At 9 = 0, vr the 
curvature vanishes. Indeed, as — > we get the sphere geometry, once that K{R = 0) = r~^, 
while for i? — oo we obtain a curvatureless infinite surface, like an infinite cylinder or an annulus. 
These points will be important later, whenever discussing upon the behavior of solitonic and 
vortex energies at these limits. Although we shall treat explicitly the case of a ring torus, 
< r < R, our approach may be extended to the other standard tori, unless otherwise specified. 






Figure 2: Global and cross-section views of a horn (upper) and spindle tori (in this case, the torus was 
cut for improving global view). The parameters are those from Fig. [TJ Note that in the cross-sections 
the horn and the spindle tori have a single and a pair of self-intercepting points, respectively. Indeed, 
in the case of the spindle torus such points correspond to two circumferences along which this surface 
crosses itself. 



In this set of coordinates the Hamiltonian ([21) gets the form below {d^ = and dg = ^): 



H = J 



-IT Jo 



R + r sin 9 



[(1 + A sin2 e) {d^Qf + sin2 9 [d^^f] + 
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+ ^^^'^''^ [(1 + A sin' e) {deef + sin' 6 (9^$)'] , (7) 

from what there follow the (static) Euler-Lagrange equations for and $, respectively: 

sin e cos ei ^ / . ^ [A {d^ef + {d^^f] + ^^'''^''^ [a {deOf + (9,$)'] 
1 it + r sm j l j 

r? — l— 7^ Sin n' 

cos (1 + A sin^ e) (deQ) + —-^ [(1 + A sin^ 6)929 + 2A sin cos 6 (deB)] + 

+ — [(1 + A sin^ 6)929 + 2A sine cos 6(9^6)1 (8) 

/? -|— Sin C' T 

COS e sin' Qde^ + -— de (sin' Qde^) + - — -ds (sin' Qds^) = (9) 

r ^ R + rsmO ^ 

As expected, the general anisotropic regime of the Heisenberg model is described by non- 
linear differential equations. Suitable non-trivial solutions can be obtained provided that some 
conditions are imposed, so that special solutions may be explicitly worked out. At this point, we 
should note that equations above resemble in form those counterparts for the planar, spherical 
and pseudospherical surfaces. Indeed, whenever R + rsin9 is identified with r, 9^sin6' or gr, 
while (j) keeps its role as the azimuth-like angle, expression above exactly recover their planar, 
spherical or pseudospherical analogues |13l[Tl]. Above, r = |r| is the planar radial distance, £H is 
the sphere radius while gr accounts for the distance measured along a pseudospherical geodesic, 
say, a hyperbole. 

3 The isotropic model and solitonic solutions 

The simplest way to seek for possible solitonic solutions associated to the present model on 
a torus is considering the isotropic regime, A = 0, and writing down the Hamiltonian ([7]), and 
its associated equations flHllQl) . in a more suitable coordinate system which allows us of getting 
the sine-Gordon equation in a simpler way (further details may be found in Refs.[T2| [T7]). Such 
a new coordinate system is parametrized by: 

asinhfecosy) asinhbsimp a cos 77 ^^^^ 

cosh b — sin 77 ' cosh b — sin rj ' cosh b — sin 7] ' 

where cp and rj vary from to 27r, while the new constant parameters a and b are both real and 
positive, what impose that the above parametrization is valid only for the ring torus, R> r. In 
terms of r and R, we have: 

R 

a= ^(R + r)(R-r) , cosh6=-, (11) 

r 

which allow the interpretation of a and b as the geometrical radius and the eccentric angle, 
respectively [17]. Conversely, we have R = a/ tanhfo and r = a/ sinhfe. In addition, it is easy to 
obtain the relations (f) = and tan(?7/2) = (R + r)/{R — r) tan(6'/2). Therefore, in this set 
of coordinates, the Hamiltonian (JTj) describing toroidally symmetric solutions, O = Q{ip) and 
$ = $(//), gets the form: 



^ ' +sinhbsin'Q(d^<i>)' 



H = J df] 
'0 Jo 



sinh b 



(12) 
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while eqs. ( ISf^ read as: 



smQcosQ{d^<^f = sinh-^bdlQ and 9^$ = 0. (13) 



The latter equation has the simplest solution = tj + tiq, which after substitution in the first 
one, gives us the sine-Gordon equation (with ^ = (fsinhb): 

9|e = sine cose, (14) 

whose the simplest solution reads [121 IlH] : 

e(^) = 2arctan(e'^^'"'^'') (15) 

with energy 

Ering = Stt J tanh('7r siuh 6) = SvrJtanh ( vry ^ — M — j (16) 



which is in agreement with the saturated Bogomol'nyi inequality [19j, Esouton = SvrJlQI, once if 
we evaluate the solitonic charge, Q = J sinQdQd^, for the solution above we exactly obtain 
Qring = tanh(7r sinh 6). Figure [3] shows how the solitonic charge behaves with the torus size, R. 
Namely, note that only when R (yo (sinhfe —>■ oo) this charge equals unity. At this limit, 
the soliton agrees with its counterpart lying in an infinite cylinder and represents a complete 
mapping from the spin sphere to the target manifold (the torus), a vr-soliton, so corresponding to 
the first homotopy class of the second homotopy group of the mapping of the spin sphere to the 
(infinite) torus, say, vr2(S'^ ^^|_r^oo) = Z. However, for finite R, such a mapping is incomplete 
and no homotopy arguments can be used for classifying solution (ITSl) as a topological excita- 
tion. Indeed, in this case, we must take into account the topology of the geometrical support: 
although the genus prevents the complete mapping from the spin sphere onto the torus (so that 
the solution presents a fractional charge, \Q\ < 1), at the same time it also ensures topological 
stability, at principle, preventing the fractional soliton from decaying against the ground-state. 
In other words, now the soliton acquires a finite characteristic length, proportional to the genus 
size, which prevents its collapse and consequently its size from vanishing. Similar scenarios are 
provided by the annulus, the truncated cone [20j, and by the punctured pseudosphere [13j . 



Once the torus is the topological product of two circles, we may wonder whether another 
solitonic solution depending on the polar angle (6* or rj), Q{9), should not also appear in this 
framework. This is indeed the case. If we reconsider Hamiltonian ([7]) with the conditions A = 0, 
e = e(^) and $ = $(0), we get: 



H2 = J 

from what we obtain: 
55$ = and 



2n 



R + r sin 6 



R + r sin 6 



(17) 



R + r sin 9 



sine cos e(9<^<l>)2 = cos 9{dge) + 



R + r sin \ 2 



(5.0) • 



The first equation is readily solved giving $ = </> + {/)o, so that the latter one is simplified to: 



i? -|- r sin ^ 



sin e cos e = cos 9{de'c)) + 



R + r sin 6 



{die) . 



(19) 
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Figure 3: How the solitonic charge associated to solution (|15p behaves in the ring torus as function 
of R/r > 1. Although Q closes the unity very fast (e.g., for R/r = 2 we obtain Qr=2t ~ 0.99995; see 
inset plot) it should be stressed that only when i? — > oo we get Q = 1. 



This expression is a generalized sine-Gordon equation, with non-constant coefficients. At prin- 
ciple, equation above admits solitonic solutions (possibly with fractional charge), but a closed 
and analytical solution is still lacking [20j. A further light into this issue may be shed if we take 
i? — > above. At this limit, equation above gets the simple form: 

sin e cos e = sin 9 cos 9{deQ) + sin^ ^(9^9) , (20) 

that is easily solved by: 

Q{9) = ±e, (21) 

which is precisely the simplest solution we have for the spherical geometry [HI [21], whose energy 
reads -Esphere = SvrJ, corresponding to the fundamental solitons with charges ±1. Therefore, we 
have seen that solitonic solutions lying on a torus generally present (or are expected to) frac- 
tional charges. These charges equal unity only at specific limits: for i? — cx) we effectively get 
a soliton on an infinite cylinder, while for ^ the spherical counterpart is obtained. 



4 Vortex-like solutions on the torus 

Geometrically, a vortex with a net winding number k 7^ 0, may be viewed as a set of spins 
rotating in a closed circuit around a core (whose center is a singular point) or a topological 
obstruction, making impossible to change the configuration to a perfectly aligned state without 
tampering with spins at an arbitrary distance from the core. Thus, once a vortex-like configura- 
tion cannot be continuously deformed to the ground-state, it acquire the status of a topologically 
stable excitation. Vortices have been intensively studied for decades in a number of physical sys- 
tems, like superfiuids and superconductors, and they have been recently observed in several 
nanomagnets as remanent states of the net magnetization [1]. In addition, they have been pro- 
posed to take part (and even to be the key elements) in some magnetoelectronic mechanisms 
with potential applications to magnetic recording, processing, sensors, and so forth [22] • 
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When we model a vortex as a continuum of spins we intend to describe only its outer re- 
gion, once inside the core the analytical treatment is expected to give only an estimate of its 
energy, shedding no light about its real structure and spins arrangement, which require nu- 
meric/simulation techniques. However, as we shall see, in the toroidal topology, a natural cutoff 
for the vortex is provided by the genus whenever i? — r > 0, so that, in a ring torus the solution 
presents no core. On the other hand, in the case of a horn torus (i? = r) a singular core takes 
place, while for R < r (spindle torus) two singular points are verified. In both cases, the cores 
appear at the self-intercepting points (discussed earlier, see Fig. |2]and related text). 

To achieve our purposes we shall consider the Planar Rotator Model, whose continuum 
Hamiltonian may be obtained from ([7j) with A = — 1 and = 7r/2 (once we are dealing with 
static solutions, our results apply equally well to the XY model), so that: 



H = J 



2tt r'2n 



JO 



, -R + r sin6' 2 



d^de . (22) 



R + r sin 9 

Demanding $ to be cylindrically symmetric, $ = $(0), Hamiltonian above yields the following: 

dl^ = =^ <l>(0) = «:0 + 0o, (23) 

where k is the charge of the vortex while 0o is a constant accounting for its global profile, 
giving no contribution to its exchange energy. The charge (vorticity) is formally defined, in the 
continuum limit, as: 

K = ^ (f f V$) ■ dt, (24) 



27r Jc V / 



where the integration is evaluated along a closed path, C, around the genus of the torus. Taking 
the solution above to the Hamiltonian fl22|) . we obtain, for a ring torus {R> r): 

i^v-ring = 47rV/t^ ^^_Z__ , (25) 

from what we see that, distinctly from other cases (planar, conical, spherical or pseudospherical 
geometries), the energy of a vortex in a ring torus does not present singularities, so that no cutoff 
needs to be introduced to prevent spurious divergences. The profile of a vortex with k = +1 
and 00 = is presented in Fig. HI Note that its energy increases with r and decreases as R is 
raised. In addition, note that as i? ^ 00 then E^^j-mg vanishes, what is expected since at this 
limit we effectively deal with the spins lying along the axial direction of an infinite cylinder, so 
being parallel to each other (in the ferromagnetic case) . Thus, instead of a true vortex we have 
the ground-state configuration in this flat geometry, whose normalized energy vanishes. 



Similar analysis may be suitably employed to the remaining supports. However, in the case 
of a horn torus (i? = r > 0), Hamiltonian (l22l) diverges for the vortex solution (1231) . This is 
not surprising since for this case the genus of the torus has been shrunk to a point (but it is 
not absent) where a singular core may be developed. Now, this spurious divergence must be 
bypassed by introducing a cutoff around the vortex core size, Iq. Doing that, the vortex energy 
on a horn torus may be easily evaluated to give: 

^v-horn = 4:71 J cot(/o/2r) , (26) 

which depends only upon the relative sizes of the core to the vortex, lo/r. The singular core 
appears exactly where the horn torus is self-intercepting. This will also happen to the spindle 
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Figure 4: Global (left) and top (right) views of a vortex with k = +1 on the ring torus. The arrows 
represent the spin field, V<I>. The genus holds as a natural cutoff, preventing any core formation where 
energy could spuriously diverge. 



torus, where a pair of singularities takes place, once this surface crosses itself at two distinct 
points (in variable 9), like follows. 

If we insist in configurations like $ = $ (0) , we clearly realize that the energy density diverges 
at those two points where the spindle torus is self-intercepting, say, i? + r sin 6* = 0. Integrating 
Hamiltonian (1221) with solution (123|1 . but keeping 6 arbitrary, gives: 



intervals. For that we should introduce suitable cutoffs around those points where energy den- 
sity given by ( l22l) blows up, say, 6'smg = arcsin(— i?/r) and vr — arcsin(— i?/r), and only later we 
should evaluate egpindie suitably in order to bypass the singularities. However, this task is very 
tedious and results in a very length expression, which will be omitted here. The main feature is 
that the well-defined vortex energy clearly shows the appearance of two singular cores located 
at both ^sing, given above. Another interesting issue is that whenever i? — >■ 0, then integration of 
Hamiltonian ( |22l) for solution (|23ll gives us exactly that result obtained in the spherical casefH]. 
namely, exhibiting the pair of cores at antipodal points. In summary, if we start off by a ring 
torus {R > r) where the vortex is coreless and decrease R, we eventually get a horn torus, where 
a core is formed. Decreasing R further a spindle torus is obtained and the vortex now presents a 
pair of singular cores, at the self-intercepting points. At the limit i? — *• these cores are located 
at diametrically opposite positions and, effectively we have a vortex lying on a sphere of radius r. 

As previously discussed for the solitonic case, we could look for 'vortex-like' solutions de- 
pending on the polar angle, say, $ = $(^). For this case, Hamiltonian (122|) yields: 



while the topological charge is given by A;^ = ^(V$) ■ dig = k'/ \JK^ — r"^ (so, valid only for 
the ring torus). Therefore, the solution of equation above may be written as: 




) • (27) 
ult above over proper 





(29) 



whose energy is easily evaluated and gives: 



Ee-ring = A-k^kIJ^—— , (30) 

which blows up as the ring torus becomes infinite, R —>■ oo. 



5 Conclusions and Prospects 

We have studied the Heisenberg exchange model for classical spins defined on a toroidal 
geometry /topology. Solitonic and vortex-type configurations were described in some details, in- 
cluding their energy, profile and behavior at some limiting cases of the torus size and geometry. 

We first considered the isotropic regime. There, only for the ring torus {R > r) solitonic 
solutions were analytically obtained. These solutions appear to bear fractional charges while 
their stability could be ensured, at principle, by the non-trivial topology of this torus, provided 
by the finite size of its central hole, i? — r > 0. As i? increases such a charge is raised but equals 
the unity only for i? — cxd (in practice, a soliton in an infinite cylinder). On the other hand, as 
i? — a soliton lying on a sphere is described. 

Now, taking the XY regime we have investigated vortex-like configurations in this support. 
In a ring torus the vortex exhibits no core and its energy density is finite everywhere. Its net 
energy vanishes as R ^ oo, once the spins are now practically parallel each other (the ferro- 
magnetic ground-state). \{ R = r (horn torus) then the vortex develops a singular core at the 
self- intercepting point, while for _R < r (spindle torus) a pair of such cores appear, so that as 
i? — > they tend to be located at antipodal points and we recover the spherical case. 

An interesting problem to be investigated is a small nanoring, say, with the geometry and 
topology of a ring torus. Actually, the expression ^^nanorings" frequently appears in the Nano- 
magnetic literature, including some possibilities for actual applications [S]. Such an object is 
generally fabricated by making a sufficiently large centered hole in a thin cylindrical nanodisk. 
Although it shares the torus topology its geometry does not, once at the disk and hole borders 
curvature changes abruptly. In addition, magnetostatic energies, which are very sensitive to 
the size and geometry of the magnet, are relevant in these cases, so that their evaluation in 
smoother supports seems to be more manageable, what justifies our prospect. Furthermore, our 
results could have some relevance for other branches were topological excitations and/or toroidal 
surfaces concern. 



Acknowledgements 

The authors thank CNPq and Fapemig for partial financial support. 



10 



References 

[1] T. Shinjo, T. Okuno, R. Hassdorf, K. Shigeto, and T. Ono, Science 289, 930 (2000); 
J. Miltat and A. Thiaville, Science 298, 555 (2002); 

A. Wachowiak, J. Wiebe, M. Bode, O. Pietzsch, M. Morgenstern, and R. Wiesendanger, 
Science 298, 557 (2002). 

[2] V.L. Berezinskii, Sov. Phys. JETP 32 493 (1970); iMd 34, 610 (1972); 
J.M. Kosterlitz and D.J. Thouless, J. Phys. C6, 1181 (1973). 

[3] E. Witten, Nucl. Phys. B 156, 269 (1979); 
G. Veneziano, Nucl. Phys. B 159, 213 (1979); 

B. AUes, M. D'Eha, and A. Di Giacomo, Nucl. Phys. B 494, 281 (1997). 

[4] See, for example, F.M.H. Crompvoets, H.L. Bethlem, R.T. Jongma, and G. Meijer, Nature 
(London) 411 (2001) 174; 

S. Gupta, K.W. Murch, K.L. Moore, T.P. Purdy, and D.M. Stamper-Kurn, Phys. Rev. Lett. 
95 (2005) 143201. 

[5] For instance, see B.I. Dunlap, Phys. Rev. B 76 (1992) 1933; 

J. Liu, H. Dai, J.H. Hafner, D.T. Colbert, R.E. Smalley, S.T. Tans, and C. Dekker, Nature 
(London) 385 (1997) 780. 

[6] S. Zhao, S. Zhang, M. Xia, E. Zhang, and X. Zuo, Phys. Lett. A331 (2004) 238. 

[7] H. Watanabe, C. Manabe, T. Shigematsu, and M. Shimizu, App. Phys. Lett. 78 (2001) 
2928. 

[8] See for example, V.P. Kravchuk, D.D. Sheka, and Y.B. Gaididei, J. Magn. Mag. Mat. 310 
(2007) 116; 

P. Vavassori, A. Busato, A. Chiapatti, A. di Bona, S. Valeri, V. Metlushko, and B. Ihc, J. 
Magn. Mag. Mat. 316 (2007) e944. 

[9] I. Hasegawa, Y. Sakaniwaa and H. Shima, J. Magn. Mag. Mat. 310 (2007) 1401; Surface 
Science 601 (2007) 5232. 

[10] A review is presented in: M.M.Hingorani and M. O' Donnel, Nature Rev. Cell Mol. Biol. 1 
(2000) 22. 

[11] M. Nakahara, Geometry, Topology and Physics, Taylor & Francis, 2nd ed. (2003). 
[12] A. Saxena, R. Dandoloff, T. Lookman, Physica A 261 (1998) 13. 

[13] L.R.A. Belo, N. M. Oliveira-Neto, W.A. Moura-Melo, A.R.Pereira, E.Ercolessi, Phys. Lett. 
A 365 (2007) 463. 

[14] G.S. Milagre, W. A. Moura-Melo, Phys. Lett. A 368 (2007) 155. 

[15] S. Wolfram, Mathematica, Cambridge Univ. Press, Cambridge, (1999). 

[16] O.-Y. Zhong-can, Phys. Rev. A 41 (1990) 4517. 

[17] J. Benoit, R. Dandoloff, Phys. Lett. A 248 (1998) 439 



11 



[18] R. Rajaraman, Solitons and Instantons: An introduction to Quantum Field Theory, North- 
Holand, Amsterdam, (1984). 

[19] E. B. Bogomol'nyi, Sov. J. Nucl. Phys. 26 (1976) 449. 

[20] A. Saxenna, R. Dandoloff, Phys. Rev. B 66, (2002) 104414. 

[21] S. Villain-Guillot, R. Dandoloff, A. Saxena, A. R. Bishop. Phys. Rev. B 52 (1995) 6712. 

[22] See, for example, M. Rahm, J. Stahl, W. Wegscheider, D. Weiss, App. Phys. Lett. 85 (2004) 
1553; 

M. Rahm, J. Stahl, D. Weiss, App. Phys. Lett. 87 (2005) 182107; 

A. R. Pereira, A. R. Moura, W. A. Moura-Melo, D. F. Carneiro, S. A. Leonel, P. Z. Coura, 
J. Appl. Phys. 101 (2007) 034310. 



12 



